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Introduction  and  summary. 

2 /Pm 

Part  Ii  The  solution  of  u„„  + u + k u • 0,  with  u (or  -q-r)  on  a 

xx  yy  ^ Z- 

b.alf-line  is  obtained  by  a simpler  and  more  direct  method  than  those 
previously  found  in  the  literature.  The  method  applies  equally  well 
to  a wedge  of  arbitrary  opening.  An  idea  of  Magnus  (1)  is  used  to 
simplify  the  solution  when  the  boundary  values  can  be  expanded  in  a 
series  of  Bessel  functions. 

Part  II;  The  problem  of  diffraction  by  an  infinite  cone  with  a 
spherical  tip  is  treated  by  a method  similar  to  the  used  in  part  I. 
Both  point  source  and  plane  wave  excitations  are  considered.  When 
the  cone  opening  is  the  problem  reduces  to  reflection  from  a 
plane  with  a hemispherical  boss.  This  simpler  problem  has  been  solved 
by  TVersky  (7):  cur*  results  agree  with  his  solution. 

Part  I. 

(a)  We  first  determine  the  Green's  function  for  the  problem. 

A two-dimensional  unit  source  is  located  at  the  point  (p*  .0» ) in 

J 

front  of  a semi- infinite  sheet  (O  ■ 0). 

We  wish  to  solve  the  two-dimensional  problem 


y 


On  leave  from  Harvard  University 
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u(^,0+)  - 0 
u(f=>,2*nr-)-  o . 


u obeys  the  Sommerfeld  radiation  condition  at  infinity 

lim  fr  (4^  - iku)  - 0 
r -xx)  r 

(r  is  measured  from  an  arbitrary  fixed  origin). 

It  is  immediately  clear  that  the  angular  eigenfunctions  of  the 

homogeneous  equation  are  just  sin  ^ ©.  This  i3  a complete  orthogonal 

set  of  functions  over  0 < © < 2 TT,  satisfying  the  boundary  conditions . 

oo 

We  writs  u(j_>,©)  - an(j=>)  sin  ^ © where  the  an(p)  are  to  be 

determined. 

We  multiply  equation  (l)  by  sin  ^ © and  integrate  from  0 to  2TT, 
obtaining: 

2 ii  2 

TTa«(f)  * J a-(f)  * k\y>)  * JS’f  sin  f de 

- sin  § ©' 

(A  priori  it  is  possible  that  the  coefficients  an(p)  are  no  better  than 

0(-™)  for  lar«:o  n,  therefore  we  may  not  be  able  to  differentiate  twice 
n 

termwise  with  respect  to  ©). 
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Integrating  twice  by  parts  and  using  the  boundary  conditions  on  u 
nA 

and  sin  we  readily  obtain 


♦ ♦ fk2  - Si 


- sin  £ O' 


JT  * (k  “ an^}  “ -/>') 


Hence 


A H^'(y)  f > j>' 
2 


B JQ(y 
2 ' 


f<r 


This  satisfies  the  radiation  condition  at  infinity  and  the  finiteness 
condition  at  the  origin. 

Ve  now  require 


\(f' 4)  • sff'  ip  *>  - of  Y * 


sin  | 


y 


Recalling  the  Wronskian  relation: 


y pp(1)(y)  Jn(y)  - ) «n(y)]  - 


We  hare: 


a„(p) 


| in  f •' 


«.  vy  r>f 

2 2 

n^(y)  Jnoy)  /<»■ 

2 2 
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We  use  the  standard  abbreviation: 


2 2 


when  ^ > yo’ 


The  symbol  > means 


i 

jj>'  when yo  </>' 


Therefore  we  have  the  Green's  function: 


(2)  u(^e)  - G(j=>,*;  - ! 5Z  sin  | e sin  § e«  (ky=>>)  Jn(kj=><) 

2 2 

(Formulas  of  this  general  nature,  derived  by  the  use  of  the  idea  of  sources 
on  a Riemann  surface,  are  to  be  found  in  A.  Sommerfeld  (2)). 

Expression  (2)  may  be  rewritten 


(ja,e; ^>',0')  - | sin  n&  sin  nS'  H^(kf>>)  Jn(kp<) 


oo 


2n+l 


h(d 


* I £ sln  nr e sin  *■  W^<> 

2 2 


n*0 


It  is  easy  to  recognize  the  first  series.  Had  we  carried  out  a similar 

analysis  for  a completely  infinite  sheet,  the  Green's  function  would 

oo  /,  \ 

turn  out  to  be  i sin  n&  sin  n£'  ;(kj®>)  Jn(kyo<). 

We  have  therefore 


^ V sin  n9  sin  n$'  H^1\kyo>)  Jn(kf=><) 


where  R is  the  distance  between  the  source  and  the  observation  point. 


♦ * 1/ 
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* 

R is  the  distance  between  the  reflected  source  about  £ * 0 and  the 
observation  point.  This  result  could  have  been  obtained  from  the 
familiar  addition  theorem  for  cylinder  functions.  We  may  write 


(3)  G (jofrjo',*')  - |^3)(kR)  - H^kR*)] 


♦ 4 y9’  sin  e sin  h£^ 

"2 


2 *»  w “ 2 v J2n+l^k/><^ 


\ 


We  are  now  in  a position  to  solve  the  problem  of  diffraction  from  a half- 

xloc 

plane.  Let  the  incoming  field  be  e propagating  towards  increasing  x. 

The  half-plane  is  located  in  the 
plane  x * 0,  with  y 0.  We  use 
the  coordinate  system  indicated 
in  Fig.  1.  The  total  field  is 

i— > x 

Fig.  1 to  vanish  for  x » 0,  y 2s  0. 

Hoc 

We  can  think  of  e as  being  the  limiting  case  of  a point  source  re- 
moved to  x * - oo. 

Indeed  suppose  there  is  a two-dimensional  unit  source  at  O’  * 


ikx 


s e 

n 


r -/*■• 


The  field  of  this  point  source 
is  just  ~ H^CkR). 


Fig.  2 


We  have: 


R 


- & sin  $)  jo1  + x 


cos 


cos 
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When  jo’  is  taken  very  large,  the  asymptotic  expansion  for 
T ^(kR)  yields 


£TT 

i F 2 “ 4 ikR 

5 JTfkR  e 6 


To  obtain  the  limit  e^*  as  j=>%  — >oo,  consider 

117  fi  a(«i 


4 

i 


"-f'  4 

r — e e 


As  <o'  -—>oo , cos — > 1,  R - yo*—> 0,  and  it  is  easy  to  see  that  the 

^ loc 

limit  of  the  above  expression  is  e 

It  is  reasonable  to  suppose  that  performing  the  same  limiting 
process  on  (3)  will  yield  the  diffraction  a plane  wave  normally  incident 
on  a half  plane.  Tbe  result  is  derived  rigorously  in  (b). 

Thus,  the  total  field  is  given  by: 


i( ft&) 


lim  ~ 

jo'-yoo 


I • 8in  nr  f 


iTT 


Tty  ^ -iy 


8 0 ■ cith*1^)*) 


H^kR*^ 


The  terms  outside  the  summation  yield 


ikx  -ikx 
e - e 


i(D 


In  the  summation, y=>  < j=>'  and  the  asymptotic  expansion  for  H^, 
is  used: 


■s&yj-J 

2 


2 JY  -in  IT 
e e 


Trk-’ 


T 
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We  obtain 


iff 


(4)  u(j=>,£)  ■ i sin  kx  + 2e 


"r  y.  (i)*n  sin  « sin  £ j 


n-0 


This  formula  can  easily  be  identified  with  that  of  Uagnus  (1)  who  uses 
the  time  dependence  e^1"^  and  a rotated  coordinate  system. 

The  results  for  the  case  where  the  normal  derivative  vanishes  on 
the  half -line  are  easily  found: 


(3a)  G(f,ej f'fi')  ♦ K^(kE*)J 


oo 


nm°  2 ~ 


.iff  oo 

(4a)  u(jo,G)  - cos  kx  + 2e  21  (i)”n  cos  i G cos  J2n+1^^ 


It  is  clear  that  the  case  of  oblique  incidence  is  easily  handled  by  the 
same  method. 

(b)  Consider  the  following  two-dimensional  boundary  value  problem: 


2 

V u + 


k2u  - 0 


u(j©,0+)  ■ u(j=>,2TT-)  ■ given  function  f^); 


u is  continuous  throughout  the  plane*  grad  u is  continuous  everywhere 
except  possibly  on  the  half-line  © ■ Oj  u obeys  a suitable  condition 
at  infinity. 


- a - 


Remark  on  the  condition  at  infinity. 

The  condition  at  infinity  is  not  the  Sommerfeld  radiation  condition 


lim  fp"  (4$  - iku)  “ 0 . 


f>cc 

In  the  case  of  diffraction  by  an  obstacle  having  some 'infinite 
dimension,  neither  the  total  field  nor  the  scattered  field  can  be  expected 
to  satisfy  the  radiation  condition. 

As  a simple  example,  consider  total  reflection  of  a plane  wave  by  an 
infinite  sheet. 


The  incident  field  is  e 
The  total  field  is 


ikx 


ikx  -ikx 
e - e 


x < 0 
x > 0 


The  "scattered"  field  is  - e“L^,XL  None  of  these  fields  obeys  the 
Sommerfeld  condition.  The  same  is  true  for  diffraction  by  a semi-infinite 
sheet. 

This  is  not  fatal,  however,  since  all  we  really  need  is 


(5) 


^ / {Bo1)‘k8>  f i-'-h i<r  ■ 0 


The  "scattered"  field  in  the  problem  of  total  reflection  obeys  condition  (5). 
In  general,  the  geometrical  optics  field  (i.e.,  no  diffraction)  will 


satisfy  (5)  C see  D.  S.  .Jones  (3)).  Since  the  scattered  field  is  expected  to 


differ  from  the  geometrical  wptics  field  only  by  a field  satisfying  the 
Sommerfeld  condition,  (5)  will  hold  for  the  scattered  field  of  the  half -plane 
problem.  In  any  event  we  adopt  (5)  as  the  imposed  condition  at  infinity.  It 
is  clear  that  if  (5)  holds  with  the  free  space  Green's  function  H^,  it 
would  also  hold  with  the  Green's  function  G of  part  I. 
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Let  us  apply  Green's  theorem  to  u and  the  Green's  function  G given 
by  (2).  As  the  boundary  of  the  region  of  integration  is  removed  to  in- 
finity, we  obtain  (in  virtue  of  (5))* 


Under  suitable  conditions,  this  may  be  written  as 


At  first,  it  might  appear  that  lim  u(^'  ,0' ) * 0;  this  is  not  true. 

o«-*o  > . 

We  will  show  that  the  integral  in  (6)  is  no  better  than  0(-4j),  hence  the 
coefficient  of  sin  — ~ O'  is  no  better  than  O(^) ; the  series  will  exhibit 
discontinuities  at  O'  * 0,  2TT,  and  lim  u(/o',0')  - f(p'). 

Formula  (6)  gives  the  solution  of  our  problem  in  a rather  untractable 
form.  Magnus  (1)  noted  that  integrals  of  the  type  (6)  are  considerably 
simplified  when  f(js>)  admits  a series  development  in  Bessel  functions  of 
integral  order.  Let 


/i00  J (kp) 

I = f m * Jl(l)  >\  T A ^ 

m,n  J u2aH{T  ' v2n+l m 


We  split  the  integral  into  integrals  from  0 to  jo ' and  yo«  to  oo 


In  each  integral  we  let  t - ky=> 
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■ {*«<*>  Wl) 


- J 


2n*l 
2 


<t)  jB<t)} 


m,n 


2 /2n+l\2 

a - 


ikp* 


**0 


H 


2 


{l' 


,d) 
*2n*r 
2 


(t)  H^(t)  - H^ft)  J, 

2 


2 /2n+l\2 


ioo 


kp« 


wy} 

2 


The  contribution  from  0 vanishes,  the  contribution  at  infinity  is 


iTT  , 2n+l\ 

~T  2 

2i  e 


IT 


*2  - (^i)2 


If  we  use  in  addition  tbe  formula  for  the  Wronskian  of  part  I,  we 
obtain 


_iTr 

. . a.  sJEah 


vm=n  . J (kP*) 

LWV)*f-JL£ 

*1\2  J 'i  2 /2t 


»,n  IT 


2 /2rv^i\ 

m - (-5-) 


<c  /2n*l\2 


For  the  case  f(j?)  ■ Jjy3)  (6)  reduces  to: 


iTT 


e“^~  x®-  n (2n+1)  sin  ^ e' 

,V“  , .J.,  w*1  > 


'vy  .«■ ) 


TT 


Vyl  2n  ♦ 1 


n-0  2 ,2a»lx2 

a - v ”o  ; 


sin  ^ O'  . 


- 11  - 


The  first  of  these  Fourier  series  in  has  coefficients  0(-^)  due 

XT 

to  the  rapid  decrease  with  n of  ^‘s  ser^es  represents  a 

continuous  function  which  vanishes  ^t  6*  ■ 0,  2TT.  Hence  the  limiting 
values  as  O'  -►  0 , 2TT - are  also  zero. 

The  second  series  has  coefficients  O(^)  and  may  exhibit  a discon- 
tinuity at  6*  * 0.  As  a matter  of  fact  it  can  be  shewn  by  an  elementary 
computation  that 

-7T  cos  m&  - 21 — "-1— - sin  0*  . 

»-0  -2  _ /2a*l*2 


Thus: 


%(/".«■)  - vy  > 


as  required.  We  can  write: 


, . sm  “iF  m (2n+l)  sin  e« 

Vr,’®,)  - LL—  £ ur  wr1  * 003 


2 ✓2^+1^2 
m - (-J-) 


CO 

In  general  if  f(^>)  admits  an  expansion  f (j°)  - ^ Jm(kpO  we  will 


obtain 
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It  is  now  easy  to  solve  the  problem  of  diffraction  by  a half -plane. 
For  simplicity  we  take  normal  incidence.  The  total  field  may  be 

split  up  in  an  incident  field  e and  a scattered  field  u(j®,©)  satisiying 
(5).  The  values  of  u(p,©)  on  the  half -plane  are  just  f(«=>)  ■ - 1.  Since 

CD  ' 

- 1 - - ]T  ep  J2p(k^),  we  can  use  (7),  with  - - €m  (-i)m  cos  m 
( € ■ 1,  €ffl  ■ 2,  m > 1),  thus  obtaining 


u(|D',e«) 


oo 


m*0 


€m;J2m(y  } 


cos  2m©' 


JflL 

oo  4 
e 

ir 


(i)-n  (2n+l)  sin  6'  J^Cy) 

■-  J ui=0 


cos  m 


*2 . <^)2 


The  first  sum  is  recognized  as  - cos  (ko*  sin  ©’ ) * - cos  kx* . The  last 


J 


sum  can  be  evaluated  by  an  elementary  computation  (see  Uagnus  (1),  p.  175) 


£: 

m-0 


cos  m 


- € 


m 


2 f2n*lsZ 
m “ '2  ' 


2 , /2n+l  TT  \ 

SH  ain  (“  2 


Thus  we  have  a single  summation  for  u(j=>*  ) and  we  obtain 


(8)  W(o,©)  - i sin  kx  *2e 


-i  X „\-n  . 2n*l 


X,y  (i)-n  sm  ZZi  • sin  X ^(y) 


which  agrees  with  formula  (4)  of  (a). 
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It  has  been  shown  that  (8)  can  be  transformed  into  the  Sommerfeld 
solution  (eeje  Magnus  (l),  p.  179)  and  into  the  KontDrowich-Lebedev 
solution  ((4),  p.  241). 


Part  II. 

We  consider  now  the  diffraction  by  a cone  with  a spherical  tip. 

The  problem  of  the  cone 


■# 

Source 


without  a tip  has  been 
treated  by  SoUfrey  (5). 


Fig.  3 

The  cone  is  infinite  and  the  source  is  taken  at  & ■ tt,  r • r'.  The 
problem  can  be  solved  with  only  slight  additional  difficulty  for  other 
locations  of  the  s puree. 

The  total  field  satisfies 

V^u  + k^u  * — ^ ^ ^ (the  field  is  independent 


2TTr  sin  e 


of 


f) 


The  boundary  condition  is  taken  as  u * 0 on  both  thr  sphere  and  cone. 


(9)  h -f?  (y2  If) 


2 

r sm  v 


: ik (sin  * 4$  * k2“ 


-§(r-r')S  (3) 

nTT 2 n4  H 

*11  t Oi.u  v 


When  r / r’,  the  homogeneous  equation  can  readily  be  separated  yielding 
the  equation  for  the  $ part  of  the  separable  solutions 


- 14  - 


ji-  3^  (sin  »»')  ♦ A » - 0 »(«0)  " 0 

w(TT)  finite 

X is  a separation  constant.  If  we  let  X ■ V (V+l),  x * cos  0, 

^ [(1-x2)  ^]  + V (v>  +l)w  - 0 w(-l)  - finite 

w(cos  0Q)  ■ 0 . 

The  only  solution  which  is  finite  at  -1  is  P^C-x),  the  Legendre  function 
of  the  first  kind. 

The  eigenvalues  V are  determined  from  the  equation  P^(-  cos  0^)  « 0. 
The  functions  Pj)(-  cos  0)  from  a complete  orthogonal  set  in  the  interval 
Oq  < 0 < TT , with  weight  sin  0.  The  normalization  factor 

J P^(-  cos  0)  sin  0 d0  » 

*0 

will  be  evaluated  later. 

We  return  to  equation  (9)  and  consider  for  u an  expansion  of  the 

foruk 


u - X.  a (r)  PJ-  cos  0)  . 

v y y 


We  multiply  (9)  by  r P^C-  cos  0)  sin  0 and  integrate  from  0^  to  T\’ 
We  obtain  the  differential  equation  for  a^(r): 

(w)  £ ^ y)  * u)>  ■ ^^rr1 
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The  solution  of  (10)  may  be  written  as: 


a (r)  - » Wte>)  5>*V2(kr<> 


A 


1ZZ 


where  Z^/2  is  a cylinder  function  which  vanishes  at  r ■ a.  We  have 
taken  care  of  the  boundary  conditions  at  r * a,  r ■ 00,  ?xd  of  con- 
tinuity at  r ■ r'.  We  determined  C from  the  jump  condition 
a‘(r'+)  - a‘(r«-)  - - 


A 


A 


2 IT  r«‘ 


aA(r'+)  " ° ^(krl)-fe  ^/2(kr,)  - H^i/2(to,)  ^2} 


Hi*iy2(kr’>  f k 

c^~*i^vv2 


(kr')  - 


z^i/z^'>TZ7z} 


Hence 


C — r Wronskian 
r' 


2TTr'2 


We  write 


>’1/2 


/ 


_ U (l ) 

>+l /2V~' 


_ca)  . 


It  is  known  that 


Wronskian  -^H^(kr),  H^(kr)J  - - ■ . 
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Hence: 


C 


1 1Kl/2<fa> 
8V  ' 


We  may  now  write  the  total  field: 


,(2) 


,(D 


u(r,e) 


i V P/-C0S  e)  Vl/2(ka)  vi/2(kr>> 
a l—  a . “rn 


>Wto>  & 


>1/2 


(kr<) 


The  strn  is  taken  over  all  ^ such  that  Pp(-  cos  O^)  * 0,  and  ^ 
defined  by  (3). 

We  now  consider  the  diffraction  of  a plane  wave  by  our  obstacle. 
This  is  to  be  obtained  by  removing  the  point  source  to  infinity  in  a 
suitable  manner.  Using  the  symbols  in  figure  2 (with  y=>»  » r'),  and 
letting  the  source  go  to  - oo  along  the  x axis: 


ikS  ikr*  ikx 
e ^ e e 

4tTH  4TTr' 


The  total  field  v(r,0)  from  a plane  wave  c 


ikx 


will  be: 


, 4TTr'  ' -* 

r >00  e 


lim 


4TTr'  v" 


r '-yOD  e 


ikr 1 


XL- 


P^C-cos 


e) 


H,(>l/2(ka> 
H 


UJ 

i+1/2 


Vl/2(kr) 


H 


(1) 


Jkr') 


(ka)  fr 


& 
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Asymptotically 


H^V2<kr‘) 


nt  e 


i(kr«  - 


VTT 


- ^L) 
2 2 ' 


The  total  field  is: 


,(2) 


(12) 


T<r,»)-Ep'(re) 

A»  Hf1.?  ,.(ka)  -fF 


(kr) 


TT_ 

2k 


^+1/2V 


may  be  obtained  from  (see  Sollfrey  (5)): 


A^  * J F^(-cos  O)  sin  © dO 


“2v-nsin  ^{^r-eos  e)  p/-cos  pv(-cos  e)j 


In  the  case  0 • is  determined  from  Pv.(0)  * 0.  Hence  ^ is  an 

o2*  y 

arbitrary  odd  integer,  2n+l 


»TT 


A2a*l  - PLl(-cos  sin  & d®  ‘ I PLn(x)  41 

1 1 W 

2 


4n+3 


v(r,e) 


cd  H^2)  ..(*&)  Z.  . ,_fka>  

V"  (An+3)  *>  ' ---  0)  —1  I JL  t 

3)‘=a*1' H^oca)  rr 


LVTT 

2 


0~G- 

o 


.4 


This  was  obtained  by  Twer sky  (7)  using  a combination  of  an  image  method 
and  superposition. 
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